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ABSTRACT 

Context. On its asteroseismic side, the initial run of CoRoT was partly devoted to the solar like star HD49933. The eigenmodes of 
this F dwarf have been observed with unprecedented accuracy. 

Aims. We investigate quantitatively the impact of changes in the modeling parameters like mass and composition. More importantly 
we investigate how a sophisticated physics affects the seismological picture of HD49933. We consider the effects of diffusion, rotation 
and the changes in convection efficiency. 

Methods. We use the CESAM stellar evolution code coupled to the ADIPLS adiabatic pulsation package to build secular models 
and their associated oscillation frequencies. We also exploited the hydrodynamical code STAGGER to perform surface convection 
calculations. The seismic variables used in this work are: the large frequency separation, the derivative of the surface phase shift, and 
the eigenfrequencies vv = o, n =i4 and v f=0 ,n=27- 

Results. Mass and uncertainties on the composition have much larger impacts on the seismic variables we consider than the rotation. 
The derivative of the surface phase shift is a promising variable for the determination of the helium content. The seismological 
variables of HD49933 are sensitive to the assumed solar composition and also to the presence of diffusion in the models. 
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1. Introduction 

vo : 

C*~) . The accurate measurements of oscillation frequencies by CoRoT 
C") ■ (Baglin et al. 2006) in nearby solar-like stars such as HD49933 
""xl" | (Appourchaux et al. 2008) aims at constraining their internal dy- 
[~^. ■ namics. First, the global parameters of the star (age, mass, corn- 
ed position) must be determined. For this, the seismic data comple- 
0\ , ment the information that already exists. Second, a comparison 
■ with sophisticated models is required. In this work we use an 
^ | updated version of the CESAM stellar evolution code (Morel 
• <-h , 1997). CESAM can take into account microscopic diffusion and 
1 the radiative accelerations of heavy elements. It also optionally 
J-j | includes the secular effects of rotation in the radiation zones 
. . i (Mathis & Zahn 2004). The resulting mixing slightly changes 
the temperature profile of otherwise stably stratified regions 
and modifies the structure and evolution of the star (Maeder & 
Meynet 2000, Decressin et al. 2009). 

Like helioseismology for the Sun, asteroseismology can con- 
strain the helium fraction in solar analogs. In this respect stars 
slightly more massive than the Sun are especially interesting tar- 
gets. Asteroseismology can also constrain convection: the depth 
of the outer convection zone depends on the efficiency of convec- 
tive energy transport in the subsurface layers and in turn affects 
the seismic signal. In the current work, we compute the relation 
between the effective temperature and the specific entropy in the 
deep convection zone where the motions become adiabatic. To 
this purpose, we use the STAGGER hydrodynamical code (Stein 
& Nordlund 1998). This computation gives an idea of the mag- 
nitude of the change of the mixing length main parameter a m \ t 
with the effective temperature for stars having surface conditions 
close to the solar ones like HD49933. 



Although some of our models are in fair agreement with the 
current classical and seismological observations of HD49933, 
our purpose is not to find the model that best fits them. We in- 
tend to estimate the influence of mass, composition and of differ- 
ent internal processes on the seismic indicators and the surface 
properties of HD49933. Section [2] presents our secular models 
of HD49933 and explores the impact of mass, composition, mi- 
croscopic diffusion and rotation. Section[3]presents the effects of 
a modified surface convection efficiency and conclusion follows 
in section [4] 



2. Secular models of HD49933 

We build models of HD49933 using the CESAM stellar evolu- 
tion code (Morel 1997). The models are initiated on the zero 
age main sequence. We use the NACRE compilation of nu- 
clear reaction rates (Angulo et al. 1999). The equation of state 
(OPAL2001) and high temperature opacities come from the 
OPAL group (Rogers, Swenson & Iglesias 1996, Iglesias & 
Rogers 1996). Below 5800 K we use the opacities of Ferguson 
et al. (2005). The equation of state and opacity tables were 
computed for the recent solar metal repartition advocated by 
Asplund et al. (2005). The atmosphere is buildt using the Hopf 
law (Mihalas 1978). It is connected to the envelope at the optical 
depth 20. Finally, the convection zone is fully homogeneous and 
modeled using mixing-length theory (hereafter MLT) in a for- 
malism close to Bohm-Vitense (1958) (see Piau et al. 2005 for 
the exact description of the current MLT formalism). Our Sun 
calibrated value for the main MLT parameter is ckmlt = 1-587. 

Solano et al. (2005) find [Fe/H]=-0.37. If we assume the 
Asplund et al. (2005) solar surface composition (X=0.7392, 
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Z=0.0122) and metal repartition then HD49933 should currently 
exhibit X=0.7664, Z = 5.4 1CT 3 provided the helium and met- 
als surface abundances follow the typical Galactic enrichment 
law ^ = 3 (Fernandes et al. 1998)Q The effective temperature 
found in the literature ranges from T e ff = 6780 + 130K (Bruntt 
et al. 2008) to T eff = 6700 + 65K (Gillon & Magain 2006) but 
values as low as T e g = 6600 + 130K have also been suggested 
(Lastennet et al. 2001). For a visual magnitude of m v = 5.77, 
the HIPPARCOS parallax (33.7 + 0.4mas, van Leeuwen 2007) 
and bolometric correction (BC = 0.025 + 0.005, Bessell et al. 
1998) suggest log(L/L ) = 0.53 + 0.01. The following models 
have been calibrated in luminosity in the sense that we stop the 
evolution once they reach the observed luminosity. Unless ex- 
plicitly mentioned all the models include overshooting of 0.2 
pressure scale height (hereafter H p ) that is typical of the ex- 
pected mass range of the star (Claret 2007). We investigate the 
impact of the mass, the composition and the rotation on the seis- 
mic properties of the star. The oscillation frequencies are com- 
puted using the Aarhus adiabatic pulsation package ADIPLS 
(Christensen-Dalsgaard & Berthomieu 1991). We focus our at- 
tention on three seismic quantities: the large frequency separa- 
tion (Av = Vf „ + i - Vf „) that is sensitive to the global properties 
of the star such as mass and age, the lowest degree lowest and 
highest order modes frequency currently observed (vf = o n= i4 and 
ve=o, B =n) and the derivative of the surface phase shift with fre- 
quency (/6*(v)) in order to address the helium ionization effects. 
We do not investigate the small separation or other seismic vari- 
ables. To determine the averages of these quantities, we excludes 
frequencies below a cut off of lmHz in order to be in the asymp- 
totic regime. 

2.1. Mass effects 

HD49933 being a single star, its mass estimate relies on evolu- 
tionary tracks. It is therefore sensitive not only to the approach 
chosen in the modelling (e. g. description of diffusion, rotation 
mixing, overshooting... ) but also to the input parameters of the 
models such as the metallicity and the helium fraction. For these 
reasons it is interesting to address the influence of the mass on 
the seismic signature of this star. Mosser et al. (2005) suggested 
M/M « 1.2 and we explore the effect of mass variation around 
this value. No diffusion or rotation effects are taken into account 
in the models of this subsection. The convection overshooting in 
the core is 0.2H p , however the overshooting has a very moderate 
impact on the variables we consider here. As illustrated by the 
work of Goupil et al. (2009), overshooting mostly influences the 
small seismic separation which we have not considered here. As 
an indication of the overshooting impact, we provide a 1.17 M 
model without overshooting on the last line of Table [TJ 

Table [TJ shows the values of the large separation and the fre- 
quencies of the modes of degree t-Q and radial orders n=14 
and n=27 '. These modes should respectively minimize and max- 
imize the surface effects. The recent analysis of Benomar et al. 
(2009) suggests two slightly different observational solutions 
owing to the modes identification. The first mode identifica- 
tion, named Ma by Benomar et al., corresponds to the work of 
Appourchaux et al. (2008) Av = 85.92 + 0.43juHz, v £=0 , n =i4 = 
1244.43 + 3.90yuHz and v f= o, n =27 = 2363.81 + 3.90/iHz. The 
other mode identification, Mb, provides Av = 85.81 + 0.29//Hz, 



1 It should be noted that the absolute Z/X=7.04 10~ 3 we adopt is 
smaller than the one adopted in the study of Goupil et al. (2009) 
Z/X=0.0106. Unlike us these authors consider the solar metal abun- 
dance and mixture of Grevesse & Noels (1993). 
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Fig. 1. Differences between the observed eigenfrequencies and 
the model eigenfrequencies vs. the radial order for degree modes 
( =0. The model considered is model 3, losanges and stars re- 
spectively correspond to the Ma and Mb modes identifications 
of Benomar et al. (2009). We also indicated the 2<x error bars 
from these authors. 

Table 1. Mass effect on the large frequency separation, Vf = o, n =i4 
and Vf=o, n =27. 



M/M G 


TcfF 


Av 


Vf = o,„=14 




Model 




[K] 


b"Hz] 


[A<Hz] 


[/'Hz] 




1.25 


6924 


100.1 


1480 


2780 


1 


1.2 


6746 


92.1 


1359 


2561 


2 


1.17 


6610 


86.2 


1270 


2403 


3 


1.17f 


6638 


87.1 


1280 


2424 


ibis 


Mode identification 


M A 




85.9 


1244 


2364 




M B 




85.8 


1290 


2406 





v f=0 ,n=i4 = 1290.55+6.97yuHzandv f=0 , n =27 = 2405.90+3.48//Hz 
(2cr error bars). The models in Table [TJ all have the same com- 
position: X=0.7664, Z = 5.4 10~ 3 . A 0.2 H p core convection 
overshooting is included in all models except the one mentioned 
by t where there is no overshooting. The second part of the table 
displays the current seismic observations following the two pos- 
sible mode identification Ma and Mb of Benomar et al. (2009). 
The best fit with seismic data is obtained for 1.17M . Within 
error bars, the large separation of this model agrees with iden- 
tification Ma of Benomar et al. (2009) and is marginaly above 
their identification M B . The absolute frequencies of our model 
and of this second solution agree at low and high order. Figure 
QJdisplays the difference in eigenfrequencies between our model 
and the two observational solutions for degree £=0, as a func- 
tion of the radial order. The difference is always below 40juHz. 
Note that in the solar case, predicted and observed eigenmodes 
frequencies agree at low frequencies but differ by ~ 20//HZ at 
high frequencies most likely because of surface effects (Turck- 
Chieze et al. 1997). The average difference between the model 
and the first and second observed solutions of Benomar et al. 
(2009) are respectively -32.6//Hz and 12.2jt/Hz. In both cases 
there is a clear trend in the model/observations difference with 
increasing frequency which might be due to the absence of sub- 
photospheric turbulence and/or magnetic surface effects as in the 
solar case. 
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The large separation, v^ = o, n =i4 and Vf=o A =n decreases rapidly 
with mass and with current modeling assumptions suggests a 
mass slightly below 1.2 M . The logarithmic derivative of the 
large separation with mass around 1.17M is = 2.6 which 
means a significant dependence of Av on mass (see hereafter). 
If we extrapolate this result to the average large separation ob- 
served by Benomar et al. (2009) 85.92yuHz, we find that the mass 
of HD49933 should be 1.168M . We note that the effective tem- 
perature of our 1.1 7M model is between 100 and 150 K below 
the current estimates but in agreement with the previous deter- 
mination of Lastennet et al.(2001). We recall that this model was 
obtained assuming [Fe/H]=-0.37, the solar Asplund et al. (2005) 
metal repartition and ^ = 3 as the metal vs. helium relation. In 
the next section we study the effects of composition. 



Table 2. Composition effects on the large frequency separation, 

Vf=0,n=14 and Vf=0,n=27- 



Metal 
repartition 

AZ 


T e ff 
[K] 


Av 
&iHz] 


V£=Q n=l4 

IjuHz] 


Vf=0 n=21 

&iHz] 


Model 


A05 

AK _ a 

AZ J 


6610 


86.2 


1270 


2403 


3 


GS98 

AK o 
AZ 


6230 


74.4 


1160 


2139 


4 


A05 
^ = 1 

A7 1 


6742 


90.8 


1338 


2525 


5 



2.2. Composition effects 

As recent developments in solar physics have shown, helioseis- 
mological results are very sensitive to the solar composition and 
therefore stand as key indicators of it (Turck-Chieze et al. 2004, 
Bahcall et al. 2005). On the one hand the metals are the main 
contributors to the opacities in low mass stars (Turck-Chieze et 
al. 2009, Turck-Chieze et al. 1993) which directly affects the 
structure and the seismic signal. On the other hand, because of 
the low effective temperature, the helium fraction is only acces- 
sible through seismology (Basu & Antia 1995). Yet the mass and 
the age of stars other than the Sun have to be known if one wants 
to constrain their metals or helium fractions (Basu et al. 2004). 
In the case of HD49933 we can therefore only give qualitative 
indications on the composition. 

We addressed the composition issue using the 1.1 7M mod- 
els, as the mass effects suggest this is the best fit for HD49933. 
The A05 and GS98 respectively stand for Asplund et al. (2005) 
and Grevesse & Sauval (1998) solar metal repartitions. In both 
cases the assumed metallicity is [Fe/H]=-0.37 dex but in the for- 
mer case this corresponds to a smaller absolute amount of met- 
als other than iron (mostly oxygen and carbon) than in the later 
case. If we assume the A05 metal repartition then X=0.7664 and 
Z = 5.4 10~ 3 (model 3) whereas if we assume the GS98 metal 
repartition X=0.7638, Z = 8.02 10 3 (model 4). These are the 
current and initial convection zone compositions because no mi- 
croscopic diffusion is taken into account in this subsection. As 
illustrated in Table|2]the metal repartition and helium abundance 
have significant impacts on Av although smaller than the im- 
pact of mass ( §2. 11 1. Please note that the models in Table [2] all 
have the same mass: 1.17M . The main difference between A05 
and GS98 metal repartition is the [O/Fe] ratio. The logarithmic 



derivative of the large separation with [O/Fe] is 



SInAv 



= 0.83 



<91nA[0/Fe] 

0. We mention that variations of the [O/Fe] ratio are observed in 
nearby stars, the scatter of [O/Fe] across the Galactic disk being 
the order of 0.2 dex (Edvardsson et al. 1993). 

The decrease in the global metal content from GS98 to A05 
creates a drop in opacity and therefore a somewhat higher ef- 
fective temperature at a given luminosity. This corresponds to a 
smaller radius and a higher interior temperature and is consistent 
with a larger large separation as: 



1 C R dr 
2A^~Jo 7 



2 Asplund et al. (2005) provide [0/Fe]=1.21 for the Sun and a vari- 
ation of [O/Fe] of s» -0.2 with respect to Grevesse & Sauval (1998) 
earlier work. 
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Fig. 2. Losanges: f3*(y) vs. v for the model having M* = 
1.17M and X=0.7664, Y=0.2282 and Z = 5.4 10~ 3 (model 
3). Triangles: f$*(v) vs. v for the same mass but X=0.7529, 
Y=0.2417 and Z = 5.3 10~ 3 (model 5). 
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Fig. 3./J*(v) vs. v for the model having M* = 1.2M Q , X=0.7664, 
Y=0.2282 and Z = 5.4 10~ 3 (model 2). 



Thus more metal rich models (GS98) suggest a mass higher than 
1.1 7M . Quantitatively, going from the A05 to the GS98 com- 
position lowers the large difference by ~ &fi Hz which is compa- 
rable to the mass effect between models 2 and 3 (O.O3M differ- 
ence). 

The region of the helium second ionization, located right 
below the photosphere, appears as a discontinuity to the in- 
ternal waves. As such it creates an oscillatory pattern in the 
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eigenmodes frequencies (Gough 1990) that has been used for 
some time to infer, for instance, the solar photosphere abun- 
dance in helium (Vorontsov et al. 1991). The oscillatory pat- 
tern is clearly visible in the second difference seismic variable 
SiVtji — vi,n-\ - ~2ve,n + vc,n+\ especially in stars slightly more 
massive than the Sun (Piau et al. 2005). While biyi.n is probably 
easier to determine from the observations than the surface phase 
shift a(v), this later variable and its derivatives are less sensitive 
to uncertainties in frequency determinations (Lopes et al. 1997). 

Following Lopes et al. (1997) we compute /3(v) = a(v) — Vjt- 
to extract the effect of helium in the cases of three stars: models 
2, 3 and 5 (see TablesQ]and|2]i. y6(v) is estimated for the low de- 
gree modes I = 0, 1 , 2. Then it is corrected by a linear fit in order 
to obtain /3*(v) whose mean value is 0. /?*(v) vs. v is shown on 
figures|2]and[3] The differences between models 2 and 3 account 
for the effects of mass whereas the differences between the mod- 
els 3 and 5 account for the effect of the helium fraction. Models 
2 and 3 have an helium mass fraction of Y=0.2282. Because of 
a different ^ law model 5 is richer in helium than the other 
models with Y=0.2417. Note that all our models have the same 
metallicity [Fe/H]=-0.37. 

Figure [2] illustrates the strong dependence of the peaks of 
jS*(v) and therefore /3(v) on the helium mass fraction. The in- 
crease of AY=0.0135 from model 3 to 5 increases the frequency 
of the first peak of /3*(v) above 1 mHz by « 150/vHz. The second 
peak of j6*(v) is shifted by ss 200juHz. However Figure [3] shows 
that the position of the peaks in /3*(v) is also strongly related to 
the mass. Based on the differences from model 3 to 5, the log- 
arithmic derivative of the large separation with the helium mass 
fraction is = 0.90. 

2.3. Diffusion effects 

The consequences of microscopic diffusion in the enveloppes 
of solar-like stars are clearly observed. In the case of the Sun 
it has been known for some time that it improves the agree- 
ment between theoretical and observed sound speed profiles 
(Christensen-Dalsgaard et al. 1993) and it is a necessary pro- 
cess to explain the current photospheric helium abundance (Basu 
& Antia 1995). Besides this, diffusion is the best candidate 
to explain the 7 Li abundance in Population II solar analogs 
(Piau 2008, Richard et al. 2005). Finally refined calculations 
by Turcotte et al. (1998) have shown that a significant change 
in surface composition through diffusion is expected between 
1.1M and 1.5M (the typical mass range of HD49933). We cal- 
culate the joint effects of microscopic diffusion (following the 
formalism of Burgers 1969) and radiative accelerations. We con- 
sider that the initial composition is the initial composition of our 
calibrated solar models including diffusion i.e. X=0.7195 and 
Y=0.2664 (Turck-Chieze et al. 2004). 

Because diffusion takes the metals away from the surface, 
the diffusive models for a given mass and Z/X surface ratio are 
cooler than the non diffusive models with the same Z/X surface 
ratio. For instance a 1.25M non diffusive model with Z/X = 
1.15 10~ 2 and Y = 0.2664 has an effective temperature of T e ff = 
6856K when it reaches the luminosity log(L/L Q ) = 0.53 at 1.12 
Gyr. The corresponding 1 .25M diffusive model exhibits T e ff = 
6 15 IK at the same luminosity which it reaches at 3.32 Gyr (see 
Table[3]l. Because of this effect, the diffusive models complying 
with the observed effective temperature of HD49933 are more 
massive than the non diffusive ones. Table [3] shows that they lie 
around 1.35M while the non diffusive model with the correct 
effective temperature have masses around 1 .2M . Regarding the 
seismic properties of diffusive models the same remark can be 



Table 3. Microscopic diffusion and radiative acceleration effects 
on the large frequency separation, Vf = o, n =i4 & Vf=o, n =27- 



M/M G 


T e ff 


Z/X 


Av 


Vf = o,„=14 


Vf=0,Ji=27 


Model 




[K] 




b"Hz] 


L»Hz] 


I^Hz] 




1.35 


6668 


4.1 10" 3 


94.0 


1401 


2627 


6 


1.3 


6480 


4.2 10~ 3 


85.7 


1280 


2401 


7 


1.25 


6151 


1.15 10- 2 


76.3 


1216 


2209 


8 



Note 1. The models all have the same initial composition (X=0.7195, 
Y=0.2664 and thus Z/X = 1.959 10~ 2 ) and are stopped at the same 
luminosity (logL/L Q = 0.53). Then the ages of the 1.25, 1.3 and 1.35 
M models are respectively 3.32, 1.79 and 0.54 Gyr. 



done. The diffusive model that best agrees with the observations 
is significantly more massive (1 .3M ) than the non diffusive one 
(1.17M ). Being more massive than the non diffusive models, 
the diffusive models fitting the observations are also younger: 
the 1.3M diffusive model is 1.79 Gyr whereas the 1.17M non 
diffusive model is 3.59 Gyr. 

As stated above, if we assume the Asplund et al. (2005) so- 
lar composition, then in HD49933 we have Z/X = 7.04 10 3 . 
Table [3] shows that the diffusive models that meet this condi- 
tion are too cool, or equivalently that the diffusive models in the 
right range of the effective temperature are too metal poor. The 
physical reason for this being that in the higher mass regime 
the outer convection zone becomes too shallow for convection 
to brake significantly the gravitational settling of the heavy el- 
ements. This drawback suggests that the microscopic diffusion 
(and radiative acceleration) effects are moderated by an addi- 
tional mixing process in the radiation envelope below the outer 
convection zone. There is actually a strong evidence of such 
a mixing as, with T e jf — 6780K HD49933, nearly stands to 
the middle of the so-called lithium dip (Boesgaard & Tripicco 
1986, Boesgaard & King 2002). The stars of the lithium dip 
show lithium and beryllium depletion which implies deep mix- 
ing in their radiation zones. A likely candidate to this is the shear 
turbulence induced by the angular momentum loss (Talon & 
Charbonnel 1998, Decressin et al. 2009). The occurance of such 
a process within HD49933 would brake the microscopic diffu- 
sion. Similarly the tachocline mixing (Spiegel & Zahn 1992) 
in the transition region between convection and radiation zones 
brakes the metals and helium diffusion in the Sun (Turck-Chieze 
et al. 2004). The interplay between diffusion and this process -or 
another non standard mixing phenomenon- appears necessary to 
build the observed 7 Li abundance in Population II (Piau 2008, 
Richard et al. 2005). In the next subsection we briefly address 
the effects of rotation and of the interaction between diffusion 
and rotation. 

2.4. Rotation effects 

The power spectrum of HD49933 suggests a surface rotation 
period of 3.4 days (Appourchaux et al. 2008). If, following 
Thevenin et al. (2006), we consider R* /R = 1 .43 this leads to 
an equatorial velocity V eq =21.3 kms _I . This result is in agree- 
ment with direct estimates of vsini (Solano et al. 2005, Mosser 
et al. 2005). 

We have applied the theory of Mathis & Zahn (2004) for two 
1.1 7M models evolving without angular momentum loss and 
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Table 4. Rotation and microscopic diffusion effects on the large 
frequency separation, ve=o, n =u & Vf = o, n =27- 
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T c ff 
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7.04 10- 


-3 


6613 


86.1 


1282 


2408 


20. 
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1.30 


7.04 10 


-3 


6480 


85.7 


1280 


2401 
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1.30t 


5.28 10 


-3 


6480 


85.7 


1281 


2401 
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Table 5. Surface gravity, effective temperature and deep convec- 
tion zone specific entropy for the hydrodynamical surface con- 
vection models of HD49933 and the best fit model, model 3, 
computed with CESAM. 





STAGGER 


STAGGER 


STAGGER 


CESAM 




model A 


model B 


model C 


model 3 


g[cm.s 2 ] 


2.75 10 4 


2.75 10 4 


2.70 10 4 


1.618 10 4 


T cff [K] 


5774 


5639 


5640 


6610 


Sadlerg.Kr'.g- 


'J1.9365 10 9 


1.9462 10 9 


1.9363 10 9 


2.3370 10 9 



showing equatorial velocities V eq = lOkms 1 (model 9) and 
V eq = 20kms~ 1 (model 10) at 4.35 Gyr and 3.97 Gyr respec- 
tively. These ages correspond to log(L/L Q ) = 0.53. The Mathis 
& Zahn (2004) transport equations account for the modifications 
in temperature, mean molecular weight and gravitational poten- 
tial induced by rotation. The meridional circulation is calculated 
as well. Table |4] provides information similar to the preceeding 
ones. The model 1 1 mentionned with f includes both rotation 
and microscopic diffusion. Its initial composition is X=0.7195, 
Y=0.2664(Z/X = 1.959 10" 2 ). All the other models of the Table 
have rotation but no microscopic diffusion effects and an initial 
composition X=0.7664, Z = 5.4 10~ 3 . The impact of rotation ef- 
fects on the large frequency separation is hardly visible. The log- 
arithmic derivative of the large separation with the equatorial ve- 
locity computed from models 9 and 10 is tiny ^°y V =1.1 10 -3 . 

Thus the rotation effects seem negligible. However when 
tested in diffusive models they diminish the diffusion efficiency 
and therefore help improving the issue outlined in 32.31 The 
physics used in building models 7 and 1 1 is similar except for 
rotation. Table|4]shows that they nearly are identical except their 
difference in the surface Z/X ratio. Model 1 1 is not our best fit 
to the observational constrains but it is the most sophisticated 
secular model of this work as it includes both of the significant 
diffusion and rotation effects expected in HD49933. Its large fre- 
quency separation and surface metallicity lie very close to the 
observations. Even tough model 1 1 is too cool, it points out that 
the mass of HD49933 predicted by refined models is higher than 
the mass predicted by simple models (e.g. model 3). 

3. Surface convection of HD49933 

3. 1 . Dynamical models 

We use the STAGGER code (Stein & Nordlund 1998) to inves- 
tigate the convective and the radiative energy transfer from 500 
km above the photosphere downto 2500 km below. The com- 
putational domain extends over 6000 by 6000 km horizontaly. 
The current grid has 63 points in each direction. We solve the 
fully compressible equations of hydrodynamics and use the same 
equation of state -OPAL- as in the CESAM calculations. As far 
as possible, we use the same opacities as well. Near and above 
the photosphere, the diffusion approximation for radiative en- 
ergy transport breaks down and one has to solve the radiative en- 
ergy transfer equation. We adopted the binning method to solve 
this equation (Nordlund 1982). This widely used method (see 
e.g. Ludwig et al. 2006) requires ID atmosphere model structure 
as well as monochromatic opacities. For both we use the Kurucz 
data and programs (Castelli 2005a, 2005b). Thus the opacities 
near and above the photosphere are not the OPAL opacities. The 



equation of state and opacity tables where computed for the com- 
position X=0.7664, Z = 5.4 10~ 3 (see 

We run STAGGER over the typical Kelvin-Helmholtz time- 
scale of the computational box that is 10 4 s. The effective tem- 
perature is computed using the surface radiative flux. The ef- 
ficiency of the surface convection is a priori dependent on the 
total energy flux (i.e. the T e jf) and the surface gravity. We con- 
sidered three simulations: a reference one for the Sun (model A), 
a model where the specific entropy of the deep convection zone 
was slightly increased with respect to the Sun (model B) and a 
model where the surface gravity was slightly lowered (model C). 
Model B is obtained by increasing the internal energy entering 
the lower boundary of the domain by 2 %. The gravity field in 
model B is kept to its solar surface value and constant through- 
out the simulation box. Model C is obtained by decreasing the 
surface gravity by 2 % but the specific entropy of the deep con- 
vection zone keeps its solar value of model A. 

Table [5] gives the surface conditions and the specific en- 
tropy in the deep convection zone (adiabatic regime) in the hy- 
drodynamical models A, B, C and our stellar evolution code 
best fit model. The increase of the specific entropy in the adi- 
abatic regime (s ac j) or the decrease of the surface gravity in- 
duce a decrease in the effective temperature. Around the solar 
surface conditions, the logarithmic derivative of s a( j with T e ff is 
= -0.225 between models A and B. 

If we extrapolate this to the case of HD49933 (T e ff = 
6780 K) it provides s a( jHD49933 = 1.8605 lO'erg.K-'.g- 1 . Our 
best fit model of HD49933 (model 3 see Table[TJ based on a so- 
lar calibrated value of ff m i t exhibits s ac j = 2.3370 10 erg.K _1 .g _1 . 
This suggests that the efficiency of convection increases when 
going from the Sun to the HD49933 surface effective temper- 
ature (decreasing s a( j corresponds to increasing a m \ t ). However 
the effects of changing the surface gravity are not accounted for. 

3.2. The mixing length parameter 

Our hydrodynamical models suggest that the efficiency of con- 
vection is higher in HD49933 than in the Sun. This means that 
the a m i t adopted for HD49933 should be larger. Table [6] esti- 
mates the effects of a modified ar m i t on Ay, vt=o n =i4 and ve = o n =27 
for 1.1 7M models with X=0.7664 and Z ='5.4 10 " 3 . As for 
the other models, core convection overshooting over 0.2H p is in- 
cluded. The specific entropies s a d of the deep convection zones 
in models 8, 3 and 9 are respectively: 2.2962 10 9 erg.K _1 .g -1 , 
2.3370 lO^rg.K '.g 1 and 2.4371 lO^rg.K '.g 1 . Our hydro- 
dynamical models Q3.lt qualitatively predict a typical change 
(s a dHD49933 - s a dModeiA) in s ad larger than what changes of 0.2 in 
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Table 6. Convection efficiency effect on the large frequency sep- 
aration, Vf =0 ,n=14 & Vf =0 ,„=27- 





T eff [K] 


Av[/(Hz] 


Vf=o.«=i4[A'Hz] 


Vf=0. n =27[^Hz] 


Model 


1.787 


6672 


88.9 


1316 


2484 


8 


1.587 


6610 


86.2 


1270 


2403 


3 


1.387 


6541 


83.5 


1471 


2317 


9 



Table 7. Logarithmic derivatives of the large separation with var- 
ious modeling parameters: mass, metal fraction, helium fraction, 
rotation and mixing length parameter. 



dhiAv 


dlnAv 


dklAv 


dhxAv 


fllnAv 


SlnAM 


SlnAlO/FeJ 


fllnAY 


fllnV cq 


d{na m[l 


2.6 


0.83 


0.90 


1.1 10~ 3 


0.21 



ttmit would induce. We conclude that a variation of 0.2 in a m \ t 
between the Sun and HD49933 is plausible. 

The trend in Av and frequency for an increasing efficiency 
of convection disagrees with the observations. Benomar et 
al. (2009) favored identification is Av = 85.92 + 0.43//Hz, 
vt= , n =u = 1244.43 + 3.90/zHz and v^ , n =27 = 2363.81 + 
3.90//Hz. For a model slightly below 1.2M this is consistent 
with a smaller a m i t throughout the evolution. The logarithmic 
derivative of the large separation with a m \ t computed from mod- 
els 8 and 9 is #2Av = 213. 

4. Conclusion 

In the context of the first CoRoT results, we build secular and 
dynamical models of the solar type star HD49933. We use 
the CESAM stellar evolution code, the hydrodynamical code 
STAGGER and the Aarhus ADIPLS oscillation package to com- 
pute the eigenfrequencies. The purposes of this work are to cal- 
culate the impact of mass, composition, diffusion, rotation and 
convection on classical and seismological parameters. We do 
not aim at finding the best fit to the current observational con- 
straints but to explore physical affects on the models. We focus 
on the (average) large frequency separation Av. In particular we 
give the logarithmic derivatives of Av with mass, composition, 
rotation and a m \ t in order to compare quantitavely their seismic 
influences. Table|7]sums up the corresponding results. 

We also investigate the derivative of the surface phase shift 
(J3(v)) that is sensitive to the helium content of the star. We fi- 
nally give the frequencies of the modes of orders n=14, n=27 
and degree I — as they respectively minimize and maximize 
the surface effects on oscillations. With one exception all the 
models include a core convection overshooting of 0.2 H p . In the 
case of the models with rotation the angular momentum loss is 
not taken into account. The results are as follows: 

1. Mass: Av is very sensitivity to mass (TableQ}. It is a factor ~ 
3 less sensitive to composition effects and a further factor ~ 
3 less sensitive to convection efficiency effects. The mass of 
HD49933 has to be tighly constrained if one wants to tackle 
composition or convection efficiency issues. 

2. Composition: the assumed solar composition strongly influ- 
ences the absolute metal content of HD49933 and therefore 
its seismic properties (Table [2}. If using seismology, a mass 



estimate to better than a few 0.0 1M will not be reliable as 
long as the solar composition is not settled. The derivative 
jS(v) of the surface phase shift depends on the helium content 
(Table|2]l. However the /3(v) dependence on mass shows once 
more that mass should be precisely constrained. 

3. Diffusion: diffusion effects are expected to be more impor- 
tant in HD49933 than in our Sun. The models taking dif- 
fusion into account suggest that HD49933 is more massive 
than the models that do not take diffusion into account (Table 
[3]). The diffusive models rise an issue: when having the cor- 
rect effective temperature they are too metal poor, when hav- 
ing the right metallicity they are too cool. 

4. Rotation: the rotational effects appear negligible on both 
classical and seismological parameters describing HD49933 
except for the Z/X ratio (Tabled. This point is also observed 
for the Sun. The interplay between rotation and diffusion im- 
proves the situation of purely diffusive models by braking the 
gravitational settling of heavy elements. 

5. Outer convection: our hydrodynamical calculations suggest 
that the efficiency of convection in HD49933 should be 
slightly higher than in the Sun. This contradicts the better 
agreement of small ar m i t models than large a m \ t models with 
the observed large frequency separation. The seismic obser- 
vations of HD49933 favor a decrease in the efficiency of con- 
vection with increasing mass. 

6. Our best fit model has a mass of 1.1 7M . It includes neither 
diffusion or rotation. Its initial composition is X = 0.7664 
and Z = 5.4 10~ 3 . This mass estimate is unchanged if core 
overshooting is suppressed. It is also not sensitive to the two 
current possible mode identification. 

This work points out the role of sophisticated physics in 
the seismological modeling of HD49933. It stresses the need 
for an accurate mass and composition determination before the 
observational data can be used to constrain the internal dynam- 
ics. Moreover it shows that diffusive effects and their interplay 
with rotational effects significantly change the mass and age es- 
timates of the star. Both diffusion and rotation effects have to 
be taken into account in HD49933 to draw reliable conculsions. 
Thus there is a need for sophisticated modelling of this CoRot 
target. That is especially important as with respect to its tem- 
perature the star lies in the middle of the so-called lithium dip. 
Therefore the future modelling works should include the angu- 
lar momentum loss as well as the possible rotationally induced 
shear turbulence in the upper radiation zone of the star. For the 
same reason it would be interesting to investigate the current sur- 
face abundance in lithium and beryllium from the observational 
point of view. 
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